Ballpark prediction for the hadronic light- by-light contribution to the muon (g — 2) 
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Using the momentum dependence of the dressed quark mass and the well-known formulae for the 
mass dependent quark loop contribution to the light-by-light scattering insertions, we compute the 
hadronic light-by-light contribution to the the muon anomalous magnetic moment. We ascribe for 
the first time a systematic error on the calculation. 

PACS numbers: 13.40.Em,14.60.Ef 
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The anomalous magnetic moment of the muon is one 
of the most accurately measured quantities in particle 
physics. Any deviation from its prediction in the Stan- 
dard Model of particle physics is a very promising signal 
of new physics. 

The present world average experimental value of its 
deviation from the Dirac value, i.e., a M = (gr„ — 2)/2, is 
given by a® xp = 11659208.9(6.3) x 1CT 10 [D H]. This 
impressive result is still limited by statistical errors, and 
a proposal to measure the muon (g — 2)^ to a precision 
of 1.6 x 10" 10 has recently been submitted to FNAL [3J. 

At the level of the experimental accuracy, the QED 
contributions to a M from photons and lcptons alone are 
very well known. Recently the calculation has been com- 
pleted up to the fifth order 0(a 5 em ) in the fine-structure 
constant a em , giving the result 11658471.885(4) x 10~ 10 



The main uncertainties at present in the Standard 
Model calculation for (g— 2) M originate from the hadronic 
vacuum polarization (HVP) as well as hadronic light- 
by-light scattering (HLBL) corrections. We show the 
present estimates and their uncertainties for the QED, 
HVP, HLBL, and the electroweak (EW) corrections in 
Table U 

TABLE I. Standard Model contributions to (g - 2) M . 



Contribution 


Result in 1(T 10 units 


Ref. 


QED(leptons) 


11658471.885 ±0.004 


m 


HVP (leading order) 


692.3 ±4.2 


A 


HVP (higher order) 


-9.84 ±0.07 


6 


HLBL 


11.6L4.0 


7 


EW 


15.4 ± 0.2 


m 


Total 


11659181.3 ±5.8 





The existing discrepancy between the experimental 
value for (g — 2) M and its Standard Model prediction 
stands at about 3a. 

In order to interpret the upcoming new experiment at 
FNAL, with an anticipated precision of 1.6 x 10 -10 , there 
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is an urgent need to improve on both the HVP as well as 
the HLBL contributions. The accuracy of the HVP con- 
tribution depends on the statistical error of the experi- 
mental data for the e + e~ annihilation cross-section into 
hadrons. With future experiments, in particular at BES- 
III [9], one foresees this error to quantitatively decrease. 
The HLBL cannot be directly related to any measur- 
able cross section however, and requires the knowledge 
of Quantum Chromodynamics (QCD) contributions at 
all energy scales. Since this is not known yet, one needs 
to rely on hadronic models to compute it. Such models 
introduce some systematic errors which are difficult to 
quantify. 

The main motivation of this work is to provide a ball- 
park prediction with a judicious error estimate for the 
HLBL scattering based on a duality argument between 
the hadronic degrees of freedom and the well-known 
quark loop contribution. 

Such a duality estimate can be obtained by invoking a 
particular regime of QCD where one knows how to per- 
form the quark loop integral responsible for the a^ LBL 
(Fig. []}. This is the large- N c of QCD [TQl CD] where a 
quark-hadron duality is accounted for considering that 
hadronic amplitudes are described by an infinite set of 
non-interacting and non-decaying resonances. As shown 
in Ref. [HI [T3J , the large- N c limit provides a very useful 
framework to approach this problem. 

Using the large- N c counting and also the chiral count- 
ing, it was proposed in [T^] to split the diagram of Fig. [I] 
into a set of different contributions where the numerically 
dominant contribution arises from the pseudo-scalar ex- 
change diagram shown in Fig. [2] [7] . 

The large- N c approach however has two shortcomings: 
firstly, the assumption of pion-exchange dominance im- 
plies that the remaining pieces are small enough to justify 
their omission. Although this seems reasonable [13], it 
might lead to an underestimation of the error. Secondly, 
calculations carried out in the large- N c limit demand an 
infinite set of resonances for computing any quantity. As 
such sum is not known in practice, one ends up truncating 
the spectral function in a resonance saturation scheme, 
the so-called Minimal Hadronic Approximation [13] . The 
resonance masses used in each calculation are then taken 
as the physical ones from PDG [15] instead of the cor- 
responding masses in the large- N c limit. Both problems 
might lead to large systematic errors not included so far 
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FIG. 1. Quark loop diagram with running quark mass. 




FIG. 2. Pion-pole contribution to HLBL scattering. The 
shaded blobs represent the Form Factor F„o~m~* (5?, gf)- 



[Mini- 
it was pointed out in Ref. [TC] that, in the large- N c 
framework, the Minimal Hadronic Approximation can 
be understood from the mathematical theory of Pade 
Approximants (PA) to meromorphic functions. Obey- 
ing the rules from this mathematical framework, one can 
compute the desired quantities in a model-independent 
way and even be able to ascribe a systematic error to the 
approach [18] . One interesting detail from this theory 
[19) is that given a low-energy expansion of a meromor- 
phic function, a PA sequence converges much faster than 
a rational function with the poles fixed in advance (such 
as the common hadronic models used so far for evalu- 
ating the HLBL), especially when the correct large Q 2 
behavior is imposed. 

To perform a quark-hadron duality estimate for the 
HLBL contribution to a M , we now discuss the direct con- 



tribution from the quark-loop diagram. The diagram 
with a light-quark running, with a mass of a few MeV, is 
only valid in the regime where quarks are not confined, 
so in the perturbative QCD regime. This is not the re- 
gion dominating the loop integral since we know that the 
momentum circulating the loop covers the low- and high- 
energy ranges. From quark models, one can use a con- 
stituent quark mass model with mass around 300 MeV 
as a rough estimate (see [2U] for a recent computation). 
This however, fails to reproduce the perturbative QCD 
result at high energies. One, then, desires a running mo- 
mentum dependent mass to perform such integration. 
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FIG. 3. Momentum dependent dressed quark mass, M(Q), 
obtained as a solution of the gap equation within the Dyson- 
Schwinger formalism. The horizontal line represents the muon 
mass. The band corresponds to an error of 10% on top of the 
chiral limit (solid red line). 



Ideally, lattice QCD calculations could be able to cal- 
culate such momentum dependent quark masses. As at 
present lattice QCD calculations arc not yet fully feasi- 
ble for physical pion masses, our proposal is to use the 
dressed-quark mass function M(Q) in the chiral limit 
computed within the Dyson-Schwinger equation (DSE) 
framework. As shown in Fig. [3j the solution of the gap 
equation for a renormalized dressed-quark propagator us- 
ing DSE provides us with the desired momentum depen- 
dent quark-mass function [5TJ[52]. We add a 10% relative 
error as an estimate for the extrapolation from the chiral 
limit to the physical light quark mass regime |23| . 

To determine the average momentum of the quark run- 
ning in the loop on Fig. [l] we employ for our duality 
argument a hadronic model. Such model yields the aver- 
age momenta entering the momentum dependent dressed 
mass M(Qi) and M(Q2), where Qi,2 refer to the momen- 
tum labels as indicated in Fig. [5] With this mass, one can 
use the well-known formulae for the a^ LBL for spin— 1/2 
fermions computed in [24-25]: 
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a BLBL (M(Q)) 



q—u,d,s 



' 16 J M(Qf 



-O 



m(q} 



log 2 



M(Qf 



(1) 



In Eq. (TTj) we display the result up to first order in 
(m^ l /M(Qy) but we include all terms up to the fourth 
order in our numerical calculations [7J. 

For our evaluation we will project the hadronic content 
of the quark loop of Fig. [I] onto the dominant hadronic 
piece of the HLBL. This is the pion-exchange contribu- 
tion considered in [27]. We will use, instead of a hadronic 
model for the transition Form Factors (FF) , a sequence of 
rational approximants |16j build up from the low-energy 
expansion of the pion FF obtained in [T7J after a fit to 
the experimental data I28f f3~l]. to minimize a model de- 
pendence. 

In practice that means to use for the FFs (blobs in 
Fig. [2]) the expression: 



F- 



P01 



.(«?,<£) «2)=« 



(2) 



where the free parameters are matched at low energies 
with the results in [TTJ: a is fixed by r w _ >77 = 7.82(24)eV 
from the PrimEx Collaboration [35] ; and b by a matching 
to the slope a v = 0.0324(22) [17]. The results for the av- 
erage momenta running in the quark loop in Fig. [T] using 
the FFs of Eq. ^ are shown in Table |n] where in paren- 
thesis we quote the errors from the input uncertainties. 

The convergence of the PA sequence to a meromorphic 
function is guaranteed by Pomerenke theorem 33 . The 
problem is to know how fast this convergence is and also 
how to ascribe a systematic error on each element of that 
sequence. For the particular case of a meromorphic func- 
tion (such as a Green's function in large- N c QCD), the 
simplest way of evaluating a systematic error is by com- 
paring the difference between two consecutive elements 
on the PA sequence [18] . 

In our approach to the FF, we evaluate the systematic 
error by computing a second element on the PA sequence 
[IB] and compare it with the result using Eq. §2§. The 
second element is: 



TABLE II. Collected results for the average momenta running 
in the quark loop in Fig. [I] its corresponding mass M(Q) in 
Fig. [3I and the a BLBL result in accordance to Eq. Q, Fig. [I] 
for both P?{Ql,Qi) and Pi(QlQ 2 2 ). 





(Qi) 


M(Q0GeV 


a HLBL x 1Q io 


P°(Qi,Ql) 


Qi 0.42(l)GeV 


0.186(3) 


12.5(3) 




Q 2 0.29(l)GeV 


0.224(2) 


9.0(1) 


PUQlQl) 


Qi 0.43(2)GeV 


0.184(6) 


12.7(7) 




Q 2 0.30(l)GeV 


0.222(3) 


9.2(2) 



1 — 2 GeV) dominates the contribution to a BLBL . To 
evaluate the error on our approximation we look for the 
maximum of the difference in the region up to 1 GeV 
between the -Pi(<7 2 ,<7 2 ) an d ^(^ij*? 2 ) as explained in 
Ref. [18]. Of course, this difference depends on the en- 
ergy, and grows as the energy increases. At 1 GeV, the 
relative difference is about 5%, and we take this error as 
a conservative estimate of the error on the whole low- 
energy region. We should add this error to the a BLBL 
results shown in Table ITU 




M(Q)//n„ 



P12 



777*7* (<7l ) fe) - 



s + bq\ 



d) (ql 



1 + bql 



>77j 



(3) 



the 



with four coefficients to be matched with T T 
slope a v , the curvature of the pion FF 6 T = 1.06(26) x 
10~ 3 LIS] and the first vector mesonresonance M p = 
0.7755(3) GeV [T3]. The results for the average momenta 
running in the quark loop in FigJTJusing the approximant 
of Eq. ([§ are shown in Table |II| where again the errors 
are due to the input ones. 

The similarity of the results obtained within both ap- 
proximants indicates that the low-energy region (up to 



FIG. 4. a^ LBL results using the expansion in Eq. M at or- 
der (ml/M{Qff (dashed blue) and at order (m^/M(Q) 2 ) 8 
(solid black) in terms of the running quark mass M(Q). The 
red point indicates the exact result for the point M(Q) = m^. 
The horizontal band shows the ballpark result from Eq. 



To quote a final number for a^ LBL implies to consider 
several sources of error: 

• Firstly, we have the error coming from the experi- 
mental inputs used to build up our approximants, 
which arises mainly from the fit to the experimental 
data on the FF. With the new forthcoming exper- 
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imental data on that FF at BES-III [5] we intuit 
lower input errors on our results. 

• Secondly, we have a 10% error due to the departure 
from the chiral limit shown in Fig. [3] This error is 
absolutely asymmetric and moves the M(Q) val- 
ues in Tabic [XT] upwards. That implies a decrease 
of about 15% on a BLBL (down to 10.6 and 7.6 in 
units of 10 -10 for Q\ and Q2 respectively). We are 
using the DSE to obtain the momentum dependent 
dressed quark masses. Once the lattice calculations 
will reach the physical quark mass values, this 10% 
error will disappear. 

• Furthermore, we also have a systematic error from 
the PA sequence used. We estimate this to be 
around 5%. That implies an error on a BLBL about 
±0.6 for Q x and ±0.5 for Q 2 in units 10~ 10 . Since 
both Pi (<7x)<zf) an d PiiliJ I2) gi ye almost the same 
results for a^ LBL , no extra error due to the diffcr- 

ence between them should be included on a BLBL . 
We should remark that the FFs employed here, al- 
though constrained by the experimental data, do 
not have the correct behavior when both photon 
virtualities (gi and 92) are very large where a be- 
havior ~ l/{ql +g|) is predicted [34H36], This fact 
does not affect our calculation since a BLBL is very 
largely dominated by the low-energy region. 

• The last source of error comes from the evalua- 



tion of a BLBL using the order (mf l /M(Q) 2 ) 8 in 
Eq. ([!]) instead of the full result. The difference 
is so smooth, Fig. |4j that no extra error should be 
considered so far. 

Our ballpark estimation lies then in the range: 



HLBL 



= [7.6(4) -T- 12.5(7)] x 10" 



-10 



(4) 



where the error in parenthesis is the combined system- 
atic and input errors as commented above and the two 
numbers represent the range due to the two momenta 
considered in our computations. 

In summary, we presented a ballpark estimate of the 
hadronic light-by-light scattering contribution to the 
(g — 2)^ based on a duality argument, and estimated the 
average momentum flowing in the quark loop diagram of 
Fig. [T] from a hadronic model. This average momentum 
then allowed us to calculate the momentum dependent 
quark mass in the quark-loop result from a^ LBL . Most 
of the recent phenomenological calculations of the a BLBL 
fall into the range obtained in this work ([T31 033 H3 133" 

Ei). 
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